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ABSTRACT 

We undertake another attempt towards seismic modelling of the most intensive stud- 
ied main sequence pulsators of the early B spectral type, v Eridani. Our analysis is 
extended by a requirement of fitting not only pulsational frequencies but also the com- 
plex amplitude of the bolometric flux variation, /, related to each mode frequency. 
This approach, called complex asteroseismology, provides a unique test of stellar pa- 
rameters, atmospheres and opacities. In particular, the concordance of the empirical 
and theoretical values of / we obtained for the high-order g mode opens a new gate in 
seismic studies of the main-sequence hybrid pulsators. The most intriguing and chal- 
lenging result is that whereas an agreement of the theoretical and empirical values of 
/ for the radial mode can be achieved only with the OPAL data, a preference for the 
OP tables is derived from the analysis of the high-order gravity mode. 
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1 INTRODUCTION 

Oscillation frequencies are the primary data used in seismic 
modelling. If oscillation spectra are rich and some regular- 
ities can be identified, then a substantial amount of infor- 
mation can be extracted about stellar interior and physics. 
This is the case for the Sun, solar-like stars and some white 
dwarfs, for which equidistances in frequencies or periods are 
observed. 

In recent years, we have entered the age of studying 
pulsating variables from space, by missions such as MOST, 
CoRoT, Kelper, which are delivering great and interest- 
ing results. But despite the tremendous number of detected 
pulsational frequencies, there is no main sequence pulsator 
showing regular structures in its oscillation spectrum. The 
best example is a J Scuti star, HD50844, for which at least 
1000 frequencies were detected from the CoRoT photometry 
and no clear regularity in frequency spacing is seen (Poretti 
et al. 2009). Moreover, all satellite missions observe in one 
wide photometric passband, therefore nothing can be said 
about mode identification and a support from ground-based 
observations is unavoidable. A hope for asterseismology of 
main-sequence stars can be put in the forthcoming BRITE 
mission which will make a two-passband photometry sup- 
ported by ground-based spectroscopy. 
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Therefore, for the time being, an improvement in seis- 
mic modelling of main-sequence pulsators can be rather 
achieved by incorporating other seismic tools related to each 
mode frequency. One of such tools is the relative ampli- 
tude of radiative flux perturbation at the level of the pho- 
tosphere, which is called the /-parameter. Because of the 
non-adiabaticity of pulsation this quantity is complex. The 
/-parameter is determined in subphotospheric layers, mainly 
in the pulsation driving region, whereas oscillation frequen- 
cies are defined mainly by the stellar interior. Therefore 
these two seismic probes are complementary to each other 
and fitting them simultaneously, which we called complex 
asteroseismology (Daszynska-Daszkiewicz & Walczak 2009, 
hereafter DDW09), can bring much stronger constraints on 
parameters of the model and theory. Theoretical values of / 
are obtained in the framework of linear nonadiabatic compu- 
tations of stellar pulsation and their empirical counterparts 
are derived from multicolour time series photometry. 

Asteroseismic potential of the /-parameter has been al- 
ready proven for some 5 Scuti stars (Daszynska-Daszkiewicz, 
Dziembowski, Pamyatnykh, 2003, Daszynska-Daszkiewicz 
et al. 2005a) and /3 Cephei stars (Daszynska-Daszkiewicz, 
Dziembowski, Pamyatnykh 2005, hereafter, DD05). In the 
case of 5 Scuti stars, the /-parameter appeared to be 
strongly dependent on the treatment of subphotospheric 
convection. For the B-type pulsators, a strong sensitivity to 
the metal content and the opacities has been found. Thus, 
a comparison of empirical and theoretical values of / can 
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yield information on the convective transport in the cooler 
pulsators and on opacities in the hotter ones. 

In the present paper, we make one more attempt to- 
wards seismic modelling of u Eridani, which is a B-type 
hybrid pulsator, i.e., one showing both /3 Cephei type pul- 
sations and frequencies typical for Slowly Pulsating B-type 
(SPB) stars. In our study we aim at parallel fitting of pul- 
sation frequencies and the corresponding values of the com- 
plex, nonadiabatic parameter /, taking into account insta- 
bility conditions. Such studies have been done recently by 
us for the P Cep star Ophiuchi (DDW09). From fitting 
two centroid frequencies of 6 Oph and the /-parameter of 
the radial mode, we have found a strong preference for the 
OPAL opacity tables. In the case of f Eri, we have more 
useable frequencies, so that stronger constraints can be ex- 
pected. Our main seismic modelling is based on fitting three 
p-mode frequencies. Then, we single out models which fit 
the forth p-mode frequency and two low frequencies corre- 
sponding to high-order g modes. Moreover, empirical values 
of / can be determined in a wide range of frequencies, i.e., 
from about 0.6 to 8 c/d. Another important result of the 
presented paper is identification of the mode degree, £, for 
all frequencies detected in the recent photometric multi-site 
campaign (Jerzykiewicz et al. 2005, hereafter JOS) using two 
methods. 

In Section 2, we give a brief review of u Eridani. Section 
3 is devoted to identifying of the spherical harmonic degree, 
£, for all fourteen modes using two approaches. A discus- 
sion of empirical values of the intrinsic mode-amplitude is 
also included. In Section 4, we present results of our seismic 
modelling of f Eri from fitting pulsational frequencies and 
then from reconciling corresponding values of the complex, 
nonadiabatic parameter /. Conclusions end the paper. 



2 u ERI AND ITS PULSATIONAL MODES 

u Eridani (HD 29248, V = 3.92 mag) is the most exten- 
sively studied main-sequence pulsator of early B spectral 
type (B2III). The star is a slow rotator, with the veloc- 
ity of rotation of about 6 km/s derived from the rotational 
splitting of pulsational frequencies (Pamyatnykh, Handler & 
Dziembowski, 2004). During the last few years, this variable 
has attracted much interest thanks to dedicated photomet- 
ric (Handler et al. 2004, JOS) and spectroscopic (Aerts et al. 
2004) multisite campaigns. These observations allowed to 
detect new frequencies typical for /3 Cep stars as well as en- 
tirely new two peaks in the SPB frequency domain. The fre- 
quency analysis revealed 14 pulsational frequencies in pho- 
tometry (Handler et al. 2004, JOS), while 9 of them, includ- 
ing one SPB frequency, were identified also in spectroscopy 
(Aerts et al. 2004). As a result, f Eridani has became one of 
the most multiperiodic B-type main-sequence pulsator and 
the updated oscillation spectrum (shown in Fig. 1) required 
to reclassify it from the /3 Cep to f3 Cep/SPB type. Recently, 
another /3 Cep/SPB star rivals it in the number of pulsation 
frequencies. From an analysis of MOST photometry of 7 
Peg, Handler et al. (2009) found 14 independent peaks: 8 of 
the /3 Cep type and 6 of the SPB type. 

Mode identification of the u Eri frequencies from pho- 
tometric and spectroscopic observations (Ridder at el. 2004) 
and from inspection of the oscillation spectrum (JOS) showed 




6.0 6.5 7.0 
frequency [c/d] 

Figure 1. The oscillation spectrum of u Eridani. Fourteen fre- 
quencies were detected from photometry and nine of them were 
found also in spectroscopic variations, marked as red peaks. The 
frequency region 1.2 — 5.2 c/d, in which no frequencies were de- 
tected, is not shown. 



the well known radial, fundamental mode and three £ — 1 
dipole modes, gi, pi and p2, with all components. Identifica- 
tion of the remaining modes, including two SPB frequencies, 
was not successful or has not been made. 

Subsequently, several papers were devoted to seismic 
analysis of f Eri. Pamyatnykh, Handler & Dziembowski 
(2004) encountered a problem with excitation of the high- 
order g mode, = 0.4328 c/d, as well as the high-frequency 
p modes with u > 6 c/d. As a solution, these authors 
proposed increasing the iron abundance in the driving re- 
gion; this excited the high frequency modes but not the g 
mode. From an analysis of the rotational splitting of two 
dipole modes, they deduced that rotation of v Eri is not 
uniform, with its rate in the p-gradient zone about three 
times faster than in the outer layers. Pamyatnykh, Handler 
& Dziembowski (2004) managed to find models which re- 
produced the three centroid frequencies corresponding to 
the £ = 0, pi, £ = l,gi and £ = l,pi modes but did not 
succeed in fitting the highest centroid frequency, a compo- 
nent of the £ = 1, P2 triplet. This was done by Ausseloos et 
al. (2004) but for higher values of the overshooting param- 
eter, ctov ~ 0.3, effective temperatures much lower than 3a 
of the observational value or non-standard chemical com- 
position. Moreover, the problem with the £ — 1,P2 mode 
excitation for models with a standard chemical composition 
remained. Both papers relied on one source of opacity data, 
OPAL, and on the old solar chemical composition (Grevesse 
& Noels, 1993). 

DDOS undertook the first attempt to fit also the /- 
parameter for the radial mode using both OPAL and OP 
tables. They showed that, in the case of B-type pulsators, 
/ strongly depends on metallicity and the adopted opacity 
data. A comparison of empirical values of / with theoretical 
ones corresponding to four seismic models with Qov ~ 0.0 
and 0.1 showed that no agreement is possible regardless of 
which opacity tables are used. 

More recently, seismic studies were done by Dziem- 
bowski & Pamyatnykh (2008), These authors used OPAL 
and OP data and included a new solar composition (Asplund 
et al. 2004, hereafter A04). In that paper, Dziembowski & 
Pamyatnykh (2008) abandoned the idea of iron accumula- 
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Figure 2. The discriminant as a function of (., obtained from fitting tlie photometric amplitude ratios and phase differences for 
14 frequencies of v Eri detected in the Stromgren uvy light-curves. Results are shown for four stellar models with the mass of 9.5 Mq, 
metallicity, Z=0.02 and 0.015, OP and OPAL opacities, and the microturbulent velocity in the atmosphere, = 2 and 8 km/s. 



tion in the Z— bump zone to explain the g mode excitation, 
because there is no evidence of an abundance anomaly in 
the atmosphere of u Eri. 

Effects of differential rotation in the analysis of three 
1=1 triplets of u Eri were included by Suarez et al. (2009) . 

Detection of high-order g-mode pulsation in massive 
stars (M > 7Mq) creates new challenges and opportunities 
for theory of stellar pulsation. Firstly, a problem of mode ex- 
citation in this frequency range for massive models should 
be solved. This can be achieved by improving opacity com- 
putations and more precise seismic modelling which should 
include, e.g., chemical composition of a given star. Secondly, 
multi-colour photometric data on modes allow to determine 
empirical values of / in a wide frequency range and to com- 
pare them with theoretical counterparts. This can provide 
unique constraints because a dependence of the /-parameter 
on the mode frequency and degree, £, is very different for p 
and g modes. 



3 IDENTIFICATION OF THE DEGREE £ OF 
PULSATIONAL MODES 

Although several years have elapsed since the second pho- 
tometric campaign of i/ Eridani (JOS), mode identification 



employing photometric observables for all 14 detected pul- 
sational peaks has never been undertaken. 

The uvy Stromgran photometry and radial velocity 
measurements allowed us to apply two approaches to iden- 
tify the degree i of the pulsational modes. In the first 
case, we compare observational values of the amplitude ra- 
tios and phase differences between available passpands with 
their theoretical counterparts and rely on theoretical val- 
ues of the /-parameter which result from linear nonadia- 
batic computations of stellar oscillations (Cugier, Dziem- 
bowski & Pamyatnykh, 1994). In the second method we 
make use of the amplitudes and phases themselves and a 
value of / is determined in the Least Square process to- 
gether with £. For a detailed description of this method 
see Daszynska-Daszkiewicz, Dziembowski & Pamyatnykh 
(2003) and DD05. In the case of B-type pulsators, the sec- 
ond method demands radial velocity measurements to get a 
unique identification of the degree £. As for input from stel- 
lar atmospheres, we use Kurucz (2004) models and Claret's 
limb darkening coefficients (Claret 2000) . Both methods pro- 
vide a discriminant as a function of £, which we named xt 
or Xej according to whether theoretical or empirical values 
of / were used. For definitions see DDW09. 

In this paper, we assume observational values of the 
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Figure 3. Values of x'e ^ function of £ determined from fitting amplitudes and phases of the light and radial velocity variations for 
nine frequencies detected in both photometric and spectroscopic observations of u Eridani. Results are shown for five models: from the 
centre and the four edges of the error box (Eq. 1). 



effective temperature and luminosity derived by DD05, i.e. 

logTeft = 4.346 ± 0.014; logL/L© = 3.94 ± 0.16, (1) 

To determine the values of Xt for the fourteen frequencies 
detected in photometric variations of v Eri, we used mod- 
els with various masses, luminosities, effective temperatures 
and metallicites within the error box, assuming the A04 
chemical mixture. As an example, in Fig. 2 we plot Xt 
a function of the mode degree, £, determined for four mod- 
els with the same mass, Al = 9.5Mq, but with two values of 
metallicity, Z = 0.020 and 0.015, two opacity tables, OPAL 
and OP, and two values of the microturbulent velocity in 
the atmosphere, = 2 and 8 km/s. 

It is well established that the dominant frequency, vi, 
is the radial fundamental mode. The three dipole modes 
with all components, (i/2, t'a, Vi), [ve, vi, 1^12) and {vs, vg, 
vio), were identified in earlier papers either from photomet- 
ric observables (Ridder et al. 2004) or from the oscillation 
spectrum itself (J05). As one can see from Fig. 2, these iden- 
tifications are confirmed by the method of amplitude ratios 
and phase differences. Identification of the degree £ for vs, 
vii and for two SPB modes, i/a and i/b, is as follows. For 
Us — 7.2009 c/d, Xt reaches low values at all £, except £ — 0, 
but there is some indication that ^ = 2. A similar result is 
obtained for fn = 6.7322 c/d, but here Xt has a minimum 
at £ = 3. We have to add in advance that whereas us matches 
£ = 2 in seismic pulsational models, the frequency vn can be 
associated only with i = 4. The same identifications for these 
two frequencies were proposed by Dziembowski & Pamyat- 
nykh (2008). Luckily, an unambiguous identification of £ was 
possible for the high-order g modes: fA = 0.433 c/d is an 
£ = 1 mode, whereas i/b = 0.614 c/d is an £ = 2 mode. For 
the latter, there is some probability of £ — 4, but it is much 
smaller because visibility of this mode is much worse. 

Now, let us discuss in more detail the values of Xt for 
z^i, the radial mode. Firstly, these values are much higher 



than 1 whichever model is considered. This is because for- 
mal errors of the photometric amplitudes and phases are 
underestimated (J05); this problem has been already dis- 
cussed by DD05. Another, more relevant feature is that a 
value of Xri'^i) decreases about twofold if the microturbu- 
lent velocity of 8 km/s is adopted instead of 2 km/s, with 
other parameters fixed. This result agrees with determina- 
tion of in the atmosphere of u Eri by Gies & Lambert 
(1992) who obtained £,t = 9.4 ± 2.6 within the LTE approx- 
imation. 

The second method of mode identification can be ap- 
plied for nine frequencies, including one SPB freqyency, 
which were found in multicolour photometry and in spec- 
troscopy. From line profile variations we derived the val- 
ues of radial velocity, defined as the first moment of the 
SiIII4553A spectral line. In Fig. 3, we show the x% depen- 
dence on £ for five models: for the centre and the four edges 
of the error box (Eq. 1). The £ identifications for first seven 
frequencies, vi — v-j, and the twelfth, 1/12, agree with the pre- 
vious ones. Now, let us deal with the low frequency mode, 
vb ~ 0.614 c/d. As can be seen from Fig. 3, the discrim- 
inant x'e reaches a minimum for £ — 2 and £ = 5. This 
approach excludes £ = 4, whereas the first method excluded 
£ = 5. Moreover, these higher £ modes suffer from strong 
disc-averaging effects. Thus, we regard the 1 = 2 identifica- 
tion as most probable. A summary of our identifications of 
£ by the two methods for all 14 pulsational frequencies of f 
Eri is given in Tab. 1. 

Other quantities extracted with the second approach 
are empirical values of mode amplitudes and complex, nona- 
diabatic parameter /. Strictly speaking, these are ampli- 
tudes defined as e = ei7"(i, 0), where e is the intrinsic mode 
amplitude and Y^{i, 0) is the spherical harmonic, dependent 
on the inclination angle, i. Thus, to obtain a value of e for 
nonradial modes one has to know i. In Tab. 2, we give exem- 
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Table 1. Most probable identification of the £ degree for the 
pulsational frequencies of v Eri from two approaches. 



frequency fc/dl 


photometry 
(theoretical f's) 


phot. + Vrad 

(empirical f's) 


= 5.7632828 


£=0 


£=0 


U2 = 5.6538767 


e=i 


£=1 


1^3 = 5.6200186 


l=\ 


£=1 


U4 = 5.6372470 


e=i 


£=1 


1/5 = 7.898200 


1=1 


1=1 


1/6 = 6.243847 


£=1 


£=1 


1/7 = 6.262917 


£=1,2, 


1=1,3 


1/8 = 7.20090 


€=?(2) 




vg = 7.91383 


£^3 




uiQ = 7.92992 


£='! 




uii = 6.73223 


£=?(3) 




ui2 = 6.22360 


€=1,2,3 


1=1,2,3 


UA = 0.432786 


£=1 




vb = 0.61440 


£=2,4 


€=2,5 



Table 2. Empirical values of e = eY^'"(i,0) and the /-parameter 
for frequencies detected in both photometry and spectroscopy. 
Stellar parameters of the central model in the error box (Eq. 1) 
were assumed (M = 9.9^0, logTefj = 4.346, log L/L© = 3.94) 
and two values of the microturbulent velocity, ^( , in stellar atmo- 
sphere were considered: 2 km/s (first lines) and 8 km/s (second 
lines). 



frequency |£| /h // it 



1/1=5.76328 0.0172(7) -9.28±0.30 0.77±0.30 2 31.4 

0.0174(5) -9.28±0.20 0.75±0.20 8 14.3 

V2 =5.65388 0.0087(3) -9.25±0.25 0.16±0.25 2 4.3 

0.0088(2) -9.26±0.19 0.15±0.19 8 2.3 

V3 =5.62002 0.0080(4) -9.04±0.30 0.02±0.30 2 5.6 

0.0081(3) -9.05±0.27 0.02±0.27 8 4.4 

VA =5.63725 0.0075(2) -9.12±0.16 0.52±0.16 2 1.5 

0.0075(2) -9.12±0.15 0.50±0.15 8 1.3 

1/5 =7.8982 0.0007(1) -11.31±0.80 4.55±0.80 2 0.3 

0.0007(1) -11.45±0.72 4.44±0.71 8 0.3 

1/6 =6.2439 0.0008(1) -9.09±1.23 -1.28±1.23 2 1.3 

0.0008(1) -9.16±1.20 -1.27±1.20 8 1.2 

1/7 =6.2629 0.0007(1) -7.09±0.38 1.29±0.38 2 0.2 

0.0007(1) -7.19±0.41 1.26±0.41 8 0.2 

i/12 =6.2236 0.0004(1) -4.04±0.75 4.54±0.75 2 0.2 

0.0004(1) -4.21±0.73 4.43±0.73 8 0.2 

VB =0.6144 0.0006(2) 9.37±4.14 16.60±4.14 2 2.1 

0.0006(2) 9.28±4.02 16.46±4.02 2 2.0 



plary values of |e| and (/ij,//) for the central model in the 
error box assuming two values of the microturbulent veloc- 
ity, ^t, in the atmosphere: 2 and 8 km/s. The last column 
contains the discriminant Xe- As one can see, there is a big 
decrease also in the value of xl; when atmospheric models 
with = 8 km/s are used. 



Although within linear approximation we cannot com- 
pare empirical values of e with their theoretical counter- 
parts, it is interesting to examine determination of these 
quantities. The radial mode has by far the largest intrinsic 
amplitude; it amounts to about 1.7 % of the stellar radius. 
Components of the € = 1, gi triplet {v2,i>3, V4) have compa- 
rable values of |e| between each others but the largest one 
is obtained for the i = l,m = mode and the smallest 
one, for the centroid mode. Of course, one has to remember 
that these quantities have to be divided by Y[^{i,Q). Con- 
sequently, if the star were seen close to the equator-on, then 
the 1= l,m = mode would have the largest intrinsic am- 
plitude. Values of |e| for the € = l,pi triplet are about ten 
times smaller than those for I = l,gi. The same is true for 
the only component of the I = 1,P2 triplet visible in pho- 
tometry and spectroscopy, i.e., the £= l,m = —1 mode. In 
the case of the SPB mode, fB = 0.614 c/d, a value of |e| is 
also comparable with those of the pi triplet. 

The second set of unknowns, i.e., the empirical values of 
the complex parameter /, can be compared with the results 
of the nonadiabatic pulsation computations. These quanti- 
ties will be included in seismic modelling of u Eri and dis- 
cussed in the next section. 



4 COMPLEX ASTREOSEISMOLOGY 

As in DDW09, our seismic modelling is performed in two 
stages. First, we aim at finding stellar models which re- 
produce observed values of pulsational frequencies. In the 
case of v Eri, we start by fitting three centroid frequencies: 
i'i{£ = 0, pi), i/4(€ = 1, gi) and pq{£ = 1, pi). Then, from this 
family of seismic models we try to single out those which fit 
also the forth p-mode frequency, i'9{£ = 1, P2), and the high- 
order g modes, i>A{i = 1) and i'b(£ = 2), assuming that they 
are axisymmetric (m = 0). 

In the next step we want to get an agreement between 
theoretical values of the /-parameter and the empirical ones 
obtained with the method used already in Sect. 3 to identify 
the degrees £. 

All computations were made using the Warsaw-New 
Jersey evolutionary code (written in its original version by 
B. Paczyhski) and the pulsational nonadiabatic code by 
Dziembowski (1977). We adopt two sources of opacity data: 
OP (Seaton 2005) and OPAL (Iglesias & Rogers, 1996), 
OPAL equation of state and the chemical mixture as deter- 
mined by A04. At ZAMS, we assumed standard hydrogen 
abundance of X = 0.7 and rotational velocity of 10 km/s. 
The effect of overshooting from the convective core was in- 
cluded according to the new formulation of Dziembowski & 
Pamyatnykh (2008). 

4.1 Fitting pulsational frequencies 

To find seismic models we survey a wide and dense grid 
of stellar parameters: M, Z, T^a, Oov In the left panel 
of Fig. 4, we put our seismic models of v Eri found with 
the OPAL and OP tables on the Z — aov plane. All these 
models reproduce three centroid frequencies: i^i(£ = 0,pi), 
j/4(€ = l,gi) and v^i^ = l,Pi)- Models to the left of the 
vertical thick line have all these three modes unstable. With 
big dots we indicated models which fit the fourth frequency. 
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Figure 4. On the left: location of the OP and OPAL seismic models which fit the three centroid frequencies of i/ Eri: ui{l = 0, pi), ^^{l = 
l,gl) and ^^{i = l,pi), on the Z — aov plane. Models to the left of the vertical tick line have all these three modes unstable. The big 
dots mark models fitting also the centroid frequency of the fourth triplet, u<j{£ = 1,P2)- Models marked with the star symbols are used 
for a comparison of the empirical and theoretical values of /. On the right: position of the OP and OPAL seismic models in the HR 
diagram and the observational error box of u Eri. Evolutionary tracks are described in the text. 




Figure 5. Values of frequencies of the OP seismic models around 
the high-order g modes, i/j^ and i/b, identified as £ = 1 and ^ = 2, 
respectively. The £ = I modes correspond to radial orders n = 
16, 17, 18 and the £ = 2 ones to n = 21, 22, 23. 



i/g{£ — 1,P2), but this mode is stable in these models. The 
other symbols (stars and circles) will be discussed later on. 
Our accuracy in fitting the frequency ui was about 10"^ 
c/d. For the other two frequencies, Vi and ve, we achieved 
an accuracy of lO-'^-lO"* c/d, depending on the model, but 
in most cases it was not worse than 10~^ c/d. The accuracy 
for the fourth frequency, ug{£ — 1,P2) was equal up to 10~* 
c/d. 

In the right panel of Fig. 4 we show location of the 
OPAL and OP seismic models in the HR diagram together 
with the error box of v Eri (see Eq. 1). These models satisfy 
the following relations: 



OP 



logL/Lo = 5.1364 log Teff - 18.3855, (2a) 



OPAL: log L/L© = 5.0956 log Tcfi - 18.2214. {2b) 

The evolutionary tracks, shown from ZAMS to TAMS, were 
computed adopting the OP-A04 opacities. The effect of the 
heavy elements abundance, Z, and the overshooting param- 
eter, Oov, on the evolutionary tracks is also presented. As one 
can see, models fitting four centroid frequencies are beyond 
the observational error box of Eri. The effective tempera- 
ture is about 2a of the observational value. A similar result 
was obtained by Ausseloos et al. (2004), but for much cooler 
temperatures or rather unacceptable chemical composition 
(X,Z). 

As can be seen from the left panel of Fig. 4, in general, 
for Oov ^ 0.25, models with larger core overshooting require 
lower Z, smaller masses and effective temperatures whereas 
for Qov -SJ 0.25, models with larger core overshooting re- 
quire higher Z, smaller masses and effective temperature. 
We have also found that in some cases for one value of Z 
there exist several seismic models with different values of 
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Table 3. Parameters of selected seismic models of u Eri (circles in the left panel of the Fig. 4). Columns from left to right are: model 
number, mass, M/Mq, age, t, effective temperature, logTeff, luminosity, log L/ Lq, metallicity, Z, overshooting parameter, Oov, radius, 
R/Rq, hydrogen abundance in the centre, Xc, and relative mass of the convective core, Mc/M. 



MODEL 


M [Mq] 


t [My] 


logTeff 


logL/L© 


Z 


Qov [Hp] 


R [Rq] 




Mc/M 


OPl 
OP2 
OPS 


9.6401 
8.5892 
8.0076 


17.S85 
2S.886 
29.092 


4.3S80 
4.3179 
4.3026 


3.898 
3.792 
3.714 


0.01800 
0.01604 
0.01615 


0.0312 
0.2297 
0.3531 


6.2701 
6.0836 
5.970 


0.24680 
0.26595 
0.27876 


0.1847 
0.1965 
0.2018 


OPALl 
OPAL2 
OPALS 


9.7706 
9.0546 
8.1709 


16.172 
19.9S5 
26.325 


4.3510 
4.3S72 
4.3155 


3.951 
3.879 
3.768 


0.01450 
0.01367 
0.01367 


0.0S40 
0.1677 
0.3436 


6.2717 
6.1503 
5.9837 


0.24691 
0.26078 
0.279S3 


0.1902 
0.1980 
0.2058 




V [c/d] 



Figure 6. The normalized instability parameter, r], plotted 
against the mode frequency with degrees £ = — 4. Two seismic 
models with about the same effective temperature were consid- 
ered: OPl and OPAL2 (cf. Tab. 3). 

core overshooting. This is connected with the avoided cross- 
ing phenomenon. Results obtained with the OP and OPAL 
opacities are quite similar but models computed with the 
OPAL data require smaller Z by about 0.0025. The OPAL 
seismic models have also larger masses and effective temper- 
atures which better match observational values of Tcff and 
L (cf. the right panel of Fig. 4). As an example, in Tab. 3 
we give the parameters for selected OP and OPAL seismic 
models marked with circles in the left panel of Fig. 4. 

The radial orders for the SPB modes, estimated from 
our seismic models, are n ~ 16 — 18 for i/a ~ 0.433 c/d 
and n ~ 21 — 23 for vb = 0.614 c/d. Let us assume that 
these modes are axisymmetric (m = 0) and check whether 
there is a seismic model amongst those presented in Fig. 4 
which reproduces also these high-order g modes. In Fig. 5, 
we show the frequency values of these orders for the OP 
seismic models as a function of the overshooting parameter 
Qov ■ In the top and bottom panel we plot frequencies around 
ua and vb, respectively. We can see that in each case there 
are two seismic models which reproduce these high-order g 
modes. The accuracy of these fits amounts to 10~* c/d. The 
corresponding radial orders for va are n = 17 and n — 18 
at ctov ~ 0.20 and 0.42, respectively. In the case of i^b, we 
get n — 22 and n = 23 at Oov ~ 0.20 and 0.39, respectively. 



We can see that the OP seismic model with Qqv ~ 0.20 fits 
simultaneously va and vb ■ The next coincidence is that this 
is the model close to the 0P2 model (see Tab. 3) which fits 
also the forth pressure mode, vg{£ = 1,P2)- A similar result 
has been obtained for the OPAL seismic models and with 
about the same values of Qov, but none of them coincides 
with the 0PAL3 model which fits 1/^(1 = 1,P2). 

Now, let us come back to the instability of the modes. 
In Fig. 6 we plot the normalized instability parameter, rj, as 
a function of frequency for two seismic models with about 
the same effective temperature but differing in all other pa- 
rameters and the opacity data (models OPl and OPAL2 in 
Tab. 3). These seismic models are marked with star symbols 
in Fig. 4. Let us remaind the reader that a mode is unstable 
if > 0. As we can see, the OP seismic models are more 
unstable; this results mainly from higher metallicity for the 
OPl model {Z = 0.018). In general, at the same Z, pulsa- 
tion computations with the OP tables give more instability 
for low frequency modes (the SPB modes) whereas with the 
OPAL data we get a wider instability in the high frequency 
range (the /3 Cep modes). This is because the Z— bump in 
the OP tables is located in deeper layers where the thermal 
time scale is longer. As we have already mentioned, the three 
highest amplitude centroid frequencies {vi, va, ua) are ex- 
cited in most seismic models (to the left of the vertical thick 
line in the left panel of Fig. 4) . In all our seismic models the 
frequency ug, corresponding to the £ = 1,P2 mode, is sta- 
ble. The highest values of rj for this mode are rj = —0.10 in 
the OP models and rj = —0.04 in the OPAL models. As for 
the high-order g modes, the frequency ua, identified as the 
£ = 1 mode, is stable in all seismic models, but ub, which is 
an £ — 2 mode, reaches instability very often. 



4.2 Constraints on opacities from the /-parameter 

Empirical values of / were obtained from the method of si- 
multaneous determination of / and £ (cf. Sect. 3). In Fig. 7 
we compare the empirical and theoretical values of / for all 
nine frequencies of v Eri found in both photometry and spec- 
troscopy. Theoretical values of / corresponding to the OPl 
and OPAL2 seismic models in Tab. 3 (marked also with star 
symbols in Fig. 4) are plotted as solid and dotted lines, re- 
spectively. Modes with the spherical harmonic degrees from 
^ = to 4 were considered. As one can see, there is quite 
a good agreement for the most pulsational frequencies in- 
cluding the SPB-type mode, vb ~ 0.614 c/d, which was 



8 J. Daszynska-Daszkiewicz and P. Walczak 



10 - 




-15 



4 5 
V [c/d] 



Figure 7. Comparison of the empirical and theoretical values of / in the whole range of pulsation frequencies observed in both photometry 
and spectroscopy of f Eri. The real and imaginary parts of / are plotted in the right and left panel, respectively. Theoretical values 
corresponding to the OPl and OPAL2 seismic models are marked with star symbols in Fig. 4 (cf. also Tab. 3). Modes with the degree £ 
from to 4 were considered. 



identified as ^ = 2. Unfortunately, the frequency va = 0.433 
c/d is not visible in spectroscopy. Around 5.6 — 5.8 c/d we 
have four points which overlap, because frequencies and 
i^2,V3,V4, are very close to each other. The discrepancy for 
the remaining frequencies is caused by inadequate accuracy 
of the photometric amplitudes and phases in the case of the 
low-amplitude modes, e.g., 1^12- 

Having empirical values of / in such a wide range of fre- 
quency is a big advantage. Particularly because the depen- 
dence of / on pulsational frequency and the mode degree, £, 
is different for p and g modes. In general, the /-parameter is 
determined by pulsational frequency and the shape of eigen- 
function as 

,3, 

where a = lu^ / (GM) is the dimensionless frequency, u 
is the angular frequency, uo — 2-kv, and G, A/, R have their 
usual meaning. Thus, for p modes the second term domi- 
nates and / is independent of I. If frequencies are low (g 
modes), the first term dominates and we have a strong de- 
pendence of the /-parameter on This is clearly visible in 
Fig. 7. We can see also that while values of / for the p modes 
change rather slowly with frequency, a rapid variation of / 
occurs for the g modes. 

The best accuracy of the values of / was obtained for 
the radial mode and the three I = l,gi modes. Their fre- 
quencies are very close to each other and so are their /- 
parameters (cf. Tab. 2). Therefore, further consideration will 
concern the /-parameter of the radial mode, v\. In Fig. 8 we 
show a comparison of the empirical and theoretical values of 
/ of the OP and OPAL seismic models. We consider models 
with effective temperatures within 3a of the observational 
value. As can be seen from Fig. 8, only with the OPAL ta- 



bles we can get an agreement of the /-values within the 
errors and this is achieved for lower values of logTcff and 
lower metallicity. The same result was obtained for the /3 
Cep star 9 Oph (DDW09). The OPAL model which fits four 
u Eri frequencies has the /-parameter marginally consistent 
with our empirical determination (the big dot). We do not 
know whether this result is a stroke of good luck or has 
some physical meaning. The effect of the microturbulent ve- 
locity in the atmosphere, ^t, on the values of / is negligible 
(cf. Tab. 2). We did not examine the effect of the hydrogen 
abundance, X, because, as has been shown by DDW09, in 
the case of early B-type pulsators, the /-parameter is only 
weakly sensitive to X. 

As we have mentioned above, it would be of particular 
interest to compare the theoretical and empirical values of 
/ also for the high-order g modes. Here we can do this for 
the frequency vb which was detected in both photometry 
and spectroscopy. In Fig. 9, we show such a comparison us- 
ing the same seismic models as those in Fig. 8. With circles 
we marked models fitting the frequency vb and with aster- 
isks, three OP models with the radial order n = 23 which 
are outside the 3o- error of the effective temperature. The 
result is quite intriguing. While in the case of the /3 Cep 
mode the empirical values of / agree with the OPAL com- 
putations, for the SPB mode, ub, an agreement for both real 
and imaginary part of / can be achieved with the OP data. 

Moreover, one can notice that in the case of both, the 
p modes and the high-order g mode, theoretical values of / 
of the OPAL and OP models agree better in the real parts 
than in the imaginary ones. 

Here, we have to discuss also a possibility that the fre- 
quency vb is a non-axisymmetric mode. As we mentioned 
in Sect. 2, the rotational velocity of v Eri is very slow and 
the effect of the azimuthal order, m, on / can be entirely 
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Figure 8. Comparison of the empirical values of / for the radial 
mode, i'i{£ = 0, pi), with theoretical ones corresponding to the 
OP and OPAL seismic models of f Eri. 
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Figure 9. Comparison of the empirical values of / for the g, 
1 = 2 mode, v^, with theoretical ones corresponding to the OP 
and OPAL seismic models of u Eri. The same radial orders, n, 
as in Fig. 5 were considered. Circles mark models which fit the 
frequency vb- Two OP models with the radial order n = 23 are 
outside the 3(7 error of the effective temperature, marked with 
asterisks. 



neglected. As for the frequency dependence of / our justifi- 
cation is as follows. Taking the stellar radius of i? ~ 6-R0, 
the rotational frequency is f^ot ~ 0.019 c/d. Given the high- 
order g modes, the Ledoux constant can be approximated 



by C„ 



Thus for ^ = 2 we get the rotational split- 



ting of ~ 0.016 c/d. In so narrow frequency range, 

i.e., 0.614 ± 0.016 c/d, both the real and imaginary parts of 
/ are nearly constant. Therefore, although the assumption 
that Vb is the m = mode can be invalid, it does not change 
above results. 



5 CONCLUSIONS 

We presented comprehensive asteroseimic modelling of the 
/3 Cep/SPB star v Eridani. We started from mode identifi- 
cation for all fourteen pulsation frequencies detected in the 
light variations from the last photometric campaign (JOS). 
For nine frequencies, seen also in the radial velocity vari- 
ation, we were able to apply the method of simultaneous 
determination of the degree i and the nonadiabatic param- 
eter / (DD05). We confirmed identifications of I from pre- 
vious determinations. Moreover, mode degrees of the SPB 
frequencies were unambiguously determined. We got i = \ 
for VA = 0.433 c/d and ^ = 2 for i^b = 0.614 c/d. 

In the next step, we looked for stellar models which fit 
the three centroid frequencies: h'i{£ = 0,pi), i'4{£ = l,gi) 
and i'(i{£ = l,Pi). Using opacities from the OP and OPAL 
tables and the A04 chemical mixture, we found a family of 
seismic models with different masses, effective temperatures, 
metallicities and core overshooting parameters. The OPAL 
seismic models agreed better with observational values of 
the effective temperature, Tcff, and luminosity, L. With two 
opacity data, we managed to find also models which fit the 
fourth centroid frequency, i'g{£ = 1,P2), but with worse ac- 
curacy in comparison with the first three ones; these models 
have effective temperatures about 2a cooler than the obser- 
vational value. 

Moreover, it turned out that some of our seismic models 
reproduce simultaneously also the high-order gravity modes, 
VA and Vb- This was the case for models with the overshoot- 
ing parameter, Oov, of about 0.2 and 0.4, independently of 
which opacity data were used. Additionally, the OP seismic 
model with Qov = 0.2 is close to the one which fits also the 
fourth p mode, vg. 

Then we extended our seismic analysis by adding the 
requirement of fitting the /-parameter. A comparison of 
the empirical and theoretical values of / for nine frequen- 
cies, which appeared in both photometry and spectroscopy, 
showed an overall agreement, especially for the modes which 
have the photometric amplitudes and phases determined 
with high accuracy. These nine frequencies are in the range 
of V £ (0.6,8) c/d, encompassing eight p modes and one 
high-order g mode. A very promising result is that the em- 
pirical value of / for the SPB- type mode, vb = 0.6144 c/d, 
agrees with theoretical computations. 

A more detailed analysis of the /-parameter for the ra- 
dial fundamental mode, vi , showed that a concordance could 
be achieved only with the OPAL opacity data. The same 
result was obtained by us for the /3 Cep variable Oph 
(DDW09). On the other hand, the same comparison for the 
high-order g mode, vb , favoured the OP data. This discrep- 
ancy calls for an explanation and indicates that something 
is still missing in the computations of opacity data. Com- 
plex asteroseismology of hybrid B-type pulsators can bring 
us closer to a solution. 
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